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Abstract. Starting from a Zimm model we study selfdiffusion in a solution of crosslinked monomers. We 
focus on the effects of the hydrodynamic interaction on the dynamics and the critical behaviour at the 
sol-gel-point. Hydrodynamic interactions cause the clusters' diffusion constant to depend not only on 
the cluster's size but also on the cluster's shape - in contrast to the Rouse model. This gives rise to a 

nontrivial scaling of the Kirkwood diffusion constant averaged over all clusters of fixed size n, D n ~ n~ b 
with b = l/d s — 1/2 given in terms of the spectral dimension d s of critical percolation clusters. The 
long-time decay of the incoherent scattering function is determined by the diffusive motion of the largest 
clusters. This implies the critical vanishing D c g ~ e a of the cluster-averaged effective diffusion constant at 
the gel point with exponent a = (3/2 — r + l/d s )/a. 

PACS. 64.60.Ht Dynamic critical phenomena - 61.25. Hq Macromolecular and polymer solutions; polymer 
melts; swelling 



1 Introduction 

The relaxation of long lived density fluctuations in gelling 
polymer solutions is dominated by the translational dif- 
fusion of individual clusters pQ. A broad distribution of 

, cluster sizes and shapes gives rise to a corresponding dis- 
tribution of relaxation times and a non-exponential decay 
of the incoherent scattering function S(q,t), which has 
been determined in several light scattering experiments 
|^P1H>| • A stretched exponential is observed in the sol phase 
S(q, t) ~ exp {— (t/t q ) x }, whereas at the critical point the 
decay is algebraic in time S(q,t) ~ t~ y . The time scale 
of the stretched exponential diverges as the distance to 

■ the critical point e tends to zero, t q ~ s~ z . The mean 
relaxation time or inverse effective diffusion constant, de- 
fined by the time integral over the incoherent scattering 
function is observed to diverge with a different exponent, 
D~ff ~ e~ a . The experimentally determined exponents 
CTfj scatter considerably and are summarized in Table ^ 
together with theoretical predictions from Rouse dynam- 
ics jniEllEl and Zimm dynamics (as derived below). 



Table 1. Comparison of critical exponents for the incoherent 
scattering function. 



exponent 


Zimm 


Rouse 


d 


m 


m 


X 


0.80 


1/2 


0.66 


0.3 - 0.8 


0.64 


y 


0.71 


0.18 


0.27 


0.2 - 0.3 


0.34 


z 


0.56 


2.22 


2.5 






a 


0.16 


1.82 


1.9 


0.5 - 1 


1.9 



From a theoretical point of view the central question is: 
How does the cluster diffusion constant scale with molecu- 
lar weight n at the gel point? One expects D n ~ rT h with 
an exponent b that depends on the underlying dynamical 
model. Several dynamical models are commonly discussed 
in the literature |12lifj] . For Rouse dynamics 6=1, inde- 
pendently of the topology of the clusters With clus- 
ter statistics according to three-dimensional bond perco- 
lation, which is generally expected to be applicable, this 
gives rise jZBHJ to a = 1.82. The Zimm model takes into 
account hydrodynamic interactions, but has so far been 
solved only for linear chains [Hj • The Zimm prediction for 
linear chains is b = 1/2, implying D n ~ 1/R n for phantom 
chains. It has been suggested [5] that this result should 
also hold for fractal structures, whose radius of gyration 
R n scales with the Hausdorff dimension df, giving rise to 
D n ~ n~ x l ds . A different line of argument starts from 
the Stokes-Einstein relation D(R) ~ (k B T)/(r}(R)R) to 
express the diffusion constant of a cluster of linear dimen- 
sion R in terms of its viscosity r](R). Using r)(R) ~ R k ^" ', 
where k is the critical exponent of the viscosity, one gets 
D(R) - l/R 1+k / u . Computer simulations [TT1TT] of dense 
melts obtain values of b within the range of 0.69 < b < 
0.99. We are not aware of estimates of D n from simulations 
that include solvent effects explicitly. However, in simula- 
tions of crosslinked particles, clusters of a given size n 
coexist with a large number of smaller clusters which may 
give rise to hydrodynamic interactions in much the same 
way as solvent molecules jll) . 

It is our intention here to solve the preaveraged Zimm 
model for arbitrarily shaped clusters. We thereby deter- 
mine the scaling behaviour D n ~ n~ b of cluster diffusion 
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constants at the gel point and work out its implications 
for the long-time decay of S(q, t). A comparison to exper- 
imental values will be given in the Conclusions. 



2 Zimm dynamics 

2.1 Hydrodynamic interactions 

We consider N point-like monomers, which are character- 
ized by their position vectors Ri, i — 1, . . . , N, in three- 
dimensional Euclidean space. Permanently formed cross- 
links constrain M randomly chosen pairs of particles 
(i e ,j e ), e — 1,...,M and give rise to a distribution of 
molecular clusters of all shapes and sizes. Crosslinks are 
modelled by Hookean springs in the potential energy 



e=l 



N 



i-.i=i 



The length I > plays the role of an inverse crosslink 
strength and physical units have been chosen such that 
k^T = 1. A given crosslink configuration Q = {i e ,je}^Li 
is specified by its N x N- connectivity matrix f and can 
be decomposed uniquely into K(Q) disjoint clusters Afc, 
which are the maximal subsets of monomers that are con- 
nected by crosslinks. 

We study the dynamics of crosslinked monomers in the 
presence of a solvent fluid, giving rise to hydrodynamic 
interactions between the monomers. Purely relaxational 
dynamics in an incompressible fluid is described by the 
equation of motion [T^ll??] 



dt 



R i (t) = -Y,H i , j (R i (t)-R j (t)) 



dU 

dR 3 {t) 



%(t) (2) 



with the mobility tensor 



kj 7 1 + (1 - 



(l + rr*). (3) 



The diagonal term in Q accounts for a frictional force 
with friction constant £ that acts when a monomer moves 
relative to the solvent. The non-diagonal term reflects 
the influence of the motion of monomer j on the sol- 
vent at the position of monomer i and is given by the 
Oseen tensor [131114) . Here t] s denotes the solvent viscos- 
ity, Si j the Kronecker symbol, r :— \r\, f := r/r, 1 is the 
three-dimensional unit matrix and the dagger indicates 
the transposition of a vector. Rouse dynamics is recov- 
ered, if the non-diagonal terms of the mobility matrix are 
neglected. The Gaussian thermal-noise forces »7j(£) in 
have zero mean and co- variance 



Vi (t)v W) = 2H id (Ri- R^Sit-t'). 



(4) 



Here S stands for the Dirac-delta function and the overbar 
indicates the Gaussian average over all realizations of rj. 



2.2 Preaveraging approximation 

The equation of motion (J2J is nonlinear due to the non- 
linear dependence of the mobility on particles' positions. 
A simple but uncontrolled approximation is the so-called 
preaveraging approximation that was first introduced by 
Kirkwood and Risemann |14| and Zimm |15| . In this ap- 
proximation the mobility matrix @ is replaced by its ex- 
pectation value (Hij) C q, which is computed with the equi- 
librium distribution, i.e. the Boltzmann weight ~ e~ u . 
Due to rotational invariance of the potential Q the aver- 
aged mobility matrix is a multiple of the identity matrix 
(Hij(Ri- J^)) eq = lH^., where 



Si,j) 



1 



1 



\Ri — Ri 



(5) 



cq 



In the computation of J5J, care has to be taken of the zero 
eigenvalues of the connectivity matrix, corresponding to 
the translation of whole clusters. To this end we regularize 
the potential by adding a confining term u> X^ili R i R i 
and let u > tend to zero subsequently. The average in © 
is conveniently performed via the Fourier representation 
of l/|r|, and the result 



1 \ 1/6 
\Ri-Rj\f e\J ttZX6 



lim ( [G(w)k< + [G(u)] jtj 



2[G(u 



n,3 



-1/2 



(G) 



involves the resolvent G(uS) := (r + uS) 1 of r . The limit 
u> I is taken by expanding the resolvent G(u>) — Eq/lo + 
Z + 0{lu) in terms of lu. Here Z := (1 - E )/r is the 
pseudoinverse of the connectivity matrix and Eq denotes 
the projector on the nullspace of r which is spanned by 
the vectors which are constant when restricted to any one 
cluster. More precisely, the matrix element [-Eo]i,j is given 
by the inverse number of monomers of the cluster if i and 
j are in the same cluster and zero otherwise (cf. section 
II. D. in 8 for details). Hence, the rhs of © vanishes 
for u> I whenever i and j belong to different clusters. 
Consequently, the preaveraged mobility matrix H cq shows 
correlations of different particles only if these particles are 
in the same cluster, in other words it is block-diagonal and 
within one block given by 



Si,j + (1 - Si t j)h (n 2 Tr/lZij) 



(7) 



For convenience we introduced IZ* 



2Z, 



and h(x) — \J xjix. The parameter k := \[§~pK ^/(6irr] s £) 
plays the role of the coupling constant of the hydrody- 
namic interaction. Formally setting k = in (JJJ yields 

j = < ' _1 ^ i J'' anc ^ ^ e Zimm model for gelation reduces 
to the Rouse model for gelation |S|I7||S]. 

It is well known that the Oseen tensor does not give 
rise to a positive-definite mobility matrix for all possible 
spatial configurations of monomers. This defect is cured 
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if the Rotne-Prager-Yamakawa tensor |lfiU17| is used in- 
stead. Again, the preaveraging procedure is done with a 
confining potential which is switched off afterwards. The 
function h is then given by 



h(x) = erf (y/x) 



1 1 - e" 



(8) 



It involves the error function "H"] and reduces to WxJtt as 
x 10. 

As a result of preaveraging we obtain the Zimm model 
for crosslinked polymers in solution 



N 



3=1 



dU 



oq 



with the co-variance of the thermal noise given by 



(9) 



(10) 



Since the connectivity matrix as well as H cq are block- 
diagonal, it follows that clusters move independently of 
each other. 

The Zimm equation """J is linear, hence it can be solved 
exactly. This is most conveniently done by introducing 
new coordinates R%{t) through the coordinate transfor- 
mation 



N 



R,{t) =: Y, [( H ° q : 



1/2 



Rj(t). 



(11) 



The resulting equation of motion for Ri (t) coincides with 
the one for a monomer in the Rouse model f6l"7l"8"| for 
crosslinked monomers, but with a formal connectivity ma- 
trix f := (H e 2) 1/2 r (H^) 1 / 2 . Its solution for (transformed) 
initial data Ri(to) is given by (cf. Section II. C. in [""]) 



of monomers is chosen independently with equal probabil- 
ity c/N, corresponding to mean- field percolation or ran- 
dom graphs. For c < c cr it = 1/2 there is no macroscopic 
cluster and almost all clusters are trees [201 • Furthermore 
all trees of size n are equally likely, (ii) Clusters are gen- 
erated according to three-dimensional continuum perco- 
lation, which is closely related to the intuitive picture of 
gelation, where monomers are more likely to be crosslinked 
when they are close to each other. Since continuum per- 
colation and lattice percolation are believed to be in the 
same universality class |21|. we employ the scaling descrip- 
tion of the latter [21], giving rise to a cluster-size distri- 
bution 

T n ~ n~ T exp{— n/n*} (14) 

for e«l and n — > oo with a typical cluster size n*(e) ~ 
e^ 1 / 17 that diverges as e — > 0. 

We denote by (A) the average of an observable A(Q) 
over all crosslink realizations Q. Partial averages 

(AU^T^/N-'Y^nAM)) (15) 
\ fc=i ' 

of A over all clusters with n sites will be of particular 
interest. The normalization r„ := ^A -1 Y^=i <W,n) rep- 
resents the average number of clusters with n sites per 
monomer. It is also known as the cluster-size distribution. 
By reordering the clusters, one gets the useful identity 



K 



fc=i 



(16) 



which is valid in the absence of an infinite cluster. 



3 Diffusion constants 



R l (t) = V Uij{t - t ) Rj(to) + / dt' U itJ (t - t') rj^t 1 ) 

J to 

(12) 



with uncorrelated random forces 



(13) 



and the time evolution matrix U(t) := exp(—3tr/£ 2 ). 
The solution of the Zimm equation @ is then obtained 
by inserting (TT""jl in (JTTJ. For a discussion of ([J"} in the case 
of dendrimers see "B"] . 



2.3 Disorder average 

To complete the description of the dynamic model we have 
to specify the statistical ensemble that determines the re- 
alizations of crosslinks in the macroscopic limit M — > oo, 
N — > oo with fixed crosslink concentration c = M/N. Two 
distributions of crosslinks will be considered, (i) Each pair 



The diffusion constant of a cluster A4 with Nk sites is 
commonly defined in terms of the long-time growth of the 
mean-squared displacement of its centre of mass R C u(t) '■= 
N k X T.teU k according to 



1 



D(A4) : = Hm - [R C u{t) - R CM (0)]' 



1 




.H eq . 





(17) 



The expression Ijl7|) for the diffusion constant has already 
been derived in 22 . Another diffusion constant has been 
introduced by Kirkwood [9111*2] 



oq 

N 2 Z— < -i,3 



It provides an upper bound to the former, 
D(M k ) < D{M k ) , 



(18) 



(19) 
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as can be shown by applying the Jensen-Peierls inequality 
to ifTTjl. 

The diffusion constants l|17() also determine the long- 
time behaviour of density fluctuations as described by the 
incoherent scattering function 



1 N 

S(q t t\g):= lim _y k e 1 «-[«*(*+*")- H *(*»)] 

to— »-oc N ^— ' 



A" 



£^ eM-q 2 tD(M k )} 



k=l 



N 



(20) 



in the stationary state to — ■> -co. Due to the blockdiago- 
nal structure of r and H cq , contributions from individual 
clusters simply add up in i|20fl. 



4 Results 

4.1 Analytical results for the Kirkwood diffusion 
constant 

The scaling behaviour of the cluster-averaged Kirkwood 
diffusion constant at the critical point 



D„ 



(D) 



(21) 



can be obtained from known results on random resistor 
networks. Starting from Eq. J7J) we identify IZij as the re- 
sistance measured between any connected pair of vertices 
(i,j) in a random resistor network |2J ■ It is obtained from 
the network of crosslinked monomers by identifying each 
Hookean spring with a resistor of unit magnitude. Note 
that this is an exact correspondence. We will infer the 
critical behaviour of D n from that of the p-th moment 
p$ : = n- 2 Y^^{^i,j)n\ c=Cct . t - In fact, we will only need 
the special case pi^ 1 ^ 2 * 1 for this purpose. Actually, mo- 
ments pn with p other than —1/2 generically occur in 
these types of network problems, see e.g. the case p = 1 
in For this reason and since it does not require more 
efforts, we treat straightaway the case with general p. In 
order to compute Pn\ we consider 



V n 2 r n p& = / d 3 x / dR R p P(R. 

n=2 J&3 J 



X 



e 



-P0i + (t-3)/ct 



(22) 



Here P(R,x)dR is the joint probability that two vertices 
of the random resistor network, whose relative position 
vector in R 3 is x, belong to the same cluster and that the 
resistance measured between them lies in the interval from 
R to R + dR [2HES1 In fact, coincides with Eq. (2.45) 
in [23]. Its asymptotic behaviour as e J, involves the 
first crossover exponent <p% of random resistor networks, 
which determines how typical resistances scale with dis- 
tance, R ~ {x^ 1 ^. The validity of requires that there 
is a divergence for e J. 0, i.e. that p<j>i — (r — 3)/er > 0. 



This is particularly the case for p = —1/2 and three- 
dimensional (and also mean-field) percolation exponents. 

Making a power-law ansatz for the n-dependence of p^\ 
we deduce from ^ and Eq. (3.9) in & that p^ ~ n 2pl , 
where 

b:=acj>i/2 = l/d s -l/2 (23) 

is expressed in terms of the spectral dimension d s |26l 
12*7] of critical percolation clusters. Hence, when using the 
preaveraged Oseen tensor in J7J), i.e. h{x) — x 1 / 2 , for the 
computation of the Kirkwood diffusion constant l|18[i . we 

get C,Dn = n _1 + K\n~ b with some constant A > 0. Thus 
D n shows a crossover from Rouse behaviour D n ~ n~ l for 
n < n(n) ~ t/(t — to Zimm behaviour 



(24) 



for asymptotically large n > n(K). Since the radius of 
gyration of Gaussian phantom clusters with n vertices 

-. / ,(G) fG) ^ 

scales as R„ ~ n ' f , where = 1/6 is the associ- 
ated fractal Hausdorff dimension [2H29II30| . Eq. (J2U im- 
plies D n ~ l/i? n in the final asymptotic regime for n, as 
suggested in j^j. 

Next we take a closer look at the Kirkwood diffu- 
sion constant D n for the special case of mean-field per- 
colation. Clearly, by inserting the appropriate mean-field 
value d s =4/3 |26U27j in (|23|l . one immediately arrives at 
b = 1/4. Alternatively, one can use the analytically known 
distribution 31 of distances on random trees to compute 
even the full asymptotic behaviour of ("R-ij) for p > — 2 
as n — > oo. To do so, we exploit that for mean-field perco- 
lation the average (•)„ is over all n n ~ 2 equally weighted 
labelled tree clusters with n vertices, see e.g. [S]. Hence, 
this average is in fact independent of the crosslink concen- 
tration c, and the resistance TZij reduces to the chemical 
distance between the vertices i and j . This gives 

n-l 



^-^ n p (n — p — 1)! 
p—i v J 



i p+1 [ dx x p+1 exp n [ dy ln(l - y) 
Jo Y Jo 



(25) 

which is independent of 1 < i ^ j < n. An asymptotic 
evaluation of the integral over x for n — > oo by Laplace's 
method then yields 



(^ J ) n n - °2P/ 2 r(l+p/2)nP/ 2 , 



(26) 



where T denotes Euler's gamma function. By setting p = 
— 1/2, we arrive again at the mean-field exponent 6=1/4 
for D n - 

4.2 Numerical results for the diffusion constant 

Now we turn to the averaged diffusion constant 



D n := (D) n \ 



(27) 
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Fig. 1. D„ at the gel point for mean field percolation and Fig. 3. D n at the gel point for three-dimensional bond perco- 
different hydrodynamic interaction strengths ft. lation. 
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K 
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Fig. 2. Critical exponents b, corresponding to a power- 
D n ~ n- b in Fig. E| 



, r. Fig. 4. Critical exponents b, corresponding to a power-law 
lawflt A.-n-MnFig.H 



fit 



of clusters of size n at the gel point. Assuming the critical 
scaling D n ~ n~ b for large cluster sizes n, the Jensen- 
Peierls inequality D n < D n , see <|19|1 , implies the inequal- 
ity 

6>6 (28) 

for the critical exponents. More detailed information on 
D n is obtained from numerical studies: As already men- 
tioned, for mean-field percolation the average is over 
all n n ~ 2 equally weighted labelled trees of size n and hence 
independent of c. Labelled trees of a given size have been 
generated via the Prufer-algorithm and handled with the 
LEDA-library .321 ■ The preaveraged mobility matrix J7J) 
is computed with h corresponding to the Rotne-Prager- 
Yamakawa tensor. The resistances Ttij in trees reduce to 
shortest graph-distances, which are computed with the Di- 
jkstra algorithm [321 ■ For each n = 1 . . . 750 the diffusion 
constant (|17|l is averaged over 100 trees, which turned out 
to yield a reliable estimate for D n . In Fig. ^ we plot D n 
as a function of n on a double-logarithmic scale for dif- 
ferent values of the hydrodynamic interaction parameter 
k. The exponent b is extracted by fitting the curves to a 
power law in the interval n S [700, 750]. Fig. [21 displays the 
exponent b for different n. The horizontal line marks the 
lower bound 6=1/4 for b. A sharp crossover is observed 
from the Rouse value 6=16 for k — to smaller values 



of 6 for non-zero k. The latter are close to and may be 
identical to the lower bound 1/4. 

For the generation of clusters according to three-di- 
mensional bond percolation we apply the Leath- Algorithm 
[33| • It generates a sequence {Ni}f =1 of connected clusters, 
in terms of which the disorder average is readily com- 
puted via (A) — \imL^ooL~ 1 J2f =1 A(Afi). This implies 

(A) n = lim^^ti^^MVEti $N t ,n for the av- 
erage over clusters of size n. The algorithm has been tested 
by comparing the numerical data of the cluster-size dis- 
tribution with the well-known scaling form of t„. Addi- 
tionally, for small values of n, we compare the number of 
clusters with known exact values [3D- For each generated 
cluster the resistances IZij are computed from the Moore- 
Penrose inverse Z of the connectivity matrix r - see below 
Eq. - and inserted into Q) with h corresponding to 
the Rotne-Prager-Yamakawa tensor. Due to high require- 
ments of memory for the computation of the inverse of H eq 
we restrict ourselves to cluster sizes n < 3000. For small 
cluster sizes we compute up to 100 diffusion constants for 
given n, for large cluster sizes this number is considerably 
smaller. The total number of generated clusters lies in the 
range of 2151 for k = 1 and 5567 for k — 0.01. The re- 
sults are shown in Figs.[31andQl In Fig.[31-Dri is plotted as 
a function of n on a double-logarithmic scale for different 
values of k. The exponent 6, extracted by fitting the curves 
in Fig. [31 to a power law in the interval n £ [500,3000], 
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is shown in Fig. 0] The horizontal line marks the lower 
bound b ss 0.25 for b, based on the value d s w 1.33 [221 
12 7| . Again one may conjecture that b — b. Like the Kirk- 
wood diffusion constant, D n also exhibits a crossover from 
Rouse to Zimm behaviour at a cluster size comparable to 
n(n). 



4.3 Incoherent scattering function 

Finally, we turn to the crosslink average S(q, t) := 
(S(q, t\Q)) of the incoherent scattering function 12Ufl in 
the sol phase. By reordering the contributions of clusters 
as in l|16(l , its long-time asymptotics is seen to be bounded 
from below by 

oo oc 

S(q, t)>J2 nT n e- D " q2f > nT n e- b " q2f (29) 

n— 1 n— 1 

according to the Jensen inequality and D n < D n . Using 
(|14|> . the sum over n on the RHS of (|29|) can be evaluated 
asymptotically for e <C 1 in terms of an integral. Up to a 
multiplicative constant, this yields (q 2 t)~ v s(t/t q ), where 
we introduced the typical relaxation time t q ~ q~ 2 e~ z , 
the critical exponents y := (t — 2)/b and z := b/a and the 
scaling function 

/>oo „ 

s(X):=X y daa 1 ^ cxp[-(a + \a- b )]. (30) 
Jo 

The scaling function s(A) is of order unity as A — > and 
decays like a stretched exponential s(A) ~ \ xy exp[— j\ x ] 
with exponent x := (1 + 6) _1 and some constant 7 > 
as A — > 00. In addition, we have verified numerically that 
the first inequality in l|29|) does not alter the long-time 

behaviour for e <C 1. Together with b = 6, as suggested by 
Figs.|2and0J this implies that the desired scaling form of 
S(q, t) for t — > 00 and e <C 1 is given by 

(q 2 t)- y s (t/t q ). (31) 

Customarily, one defines an effective diffusion constant 
-Dcff by 

D^ 1 := lima 2 / dtS(q,t) = £ nr n (1/D) fl | . 

(32) 

We conclude from the scaling form (|31Jl of S(q, t), that D ff 
vanishes at the gel point according to D c s ~ e a provided 
a := (2— r+&)/er > 0. Three-dimensional bond percolation 
leads to the value a w 0.16. If instead 6 < r — 2, then Z? G g- 
remains non-zero at the transition. Such an unphysical 
situation will occur, for example, if one chooses mean- 
field percolation for the crosslink average. Note that the 
average time as exemplified by -D e ff is not proportional 
to the time scale t q of the stretched exponential, as is 
sometimes assumed incorrectly. 



5 Conclusions 

We have computed the averaged diffusion constant D n 
of clusters of a given size n for Zimm dynamics at the 
gelation transition. Our main result, as suggested by the 

data, is the scaling D n ~ n~ b . Here the exponent b is given 
by Eq. ©, and the average is with respect to percolation 
statistics. Recalling that the radius of gyration of phantom 

clusters scales like R n ~ n b , the above relation shows that 
D„ ~ 1/R n does not only hold for linear chains [Hj but, 
in an average sense, for all percolation clusters. 

Our results pertain to ^-conditions, in so far as ex- 
cluded-volume interactions have been ignored. One exper- 
imental setup to realize percolation statistics is crosslink- 
ing of a melt close to the gel point. Subsequent dilution 
is required to observe properties of single clusters as, for 
example, D n . If, instead, one were to crosslink a dilute so- 
lution, one might generate different cluster statistics and 
hence different critical exponents. 

The above scaling of the diffusion constants determines 
the critical dynamics of density fluctuations at the gelation 
transition within the framework of the Zimm model. For 
three-dimensional bond percolation the numerical values 
of the critical exponents are summarized in Table ^ an d 
compared to the predictions of the Rouse model as well 
as to experimental values. Even though the latter scatter 
considerably, it is clear that the Zimm model fails to ac- 
count for the observed critical behaviour. The failure can 
be traced to the very slow decrease of the diffusion con- 
stant with cluster size, resulting in too weak a divergence 
of the time scale and too fast a decay of the incoherent 
scattering function. 

There are at least three reasons for the discrepancy 
between the predictions of the Zimm model and experi- 
ments: (i) It has been suggested that hydrodynamic 
interactions between monomers on a cluster are screened 
by smaller clusters in the reaction bath so that the Rouse 
rather than the Zimm model should apply. Our analy- 
sis supports this conclusion in so far as the exponents of 
the Rouse model are closer to the experimental values, 
(ii) Excluded-volume interactions cause a swelling of the 
clusters whose influence on the diffusion constants can 
be estimated from a simple scaling argument: Provided 
the relation D n ~ 1/R n still holds in the presence of 
excluded-volume interactions, the scaling R n ~ n 1 / df can 
be obtained from a standard Flory argument |SS| with 
the swollen fractal dimension df = d s (d + 2)/(d s + 2). 
Here d — 3 is the space dimension, and from d s « 4/3 
one gets 1/df w 1/2, which is roughly twice as large as 
b. The resulting exponent for the vanishing of the effec- 
tive diffusion constant D c h is a 0.71 for the case of 
three-dimensional bond percolation, respectively a = 
for mean-field percolation. Yet, these values for a still lie 
well below the experimental ones. This may be due to the 
neglect of excluded-volume interactions between different 
clusters in this simple argument, (iii) Preaveraging of the 
hydrodynamic interactions is an uncontrolled approxima- 
tion, and it remains to be seen what a full treatment of hy- 



M. Kiintzel et al.: Diffusion of gelation clusters in the Zimm model 



7 



drodynamic interactions predicts for the critical dynamics 
of gelling solutions. 
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